NON-NEGATIVE GLOBAL WEAK SOLUTIONS FOR A 
DEGENERATE PARABOLIC SYSTEM MODELING THIN FILMS 

DRIVEN BY CAPILLARITY 



BOGDAN-VASILE MATIOC 



Abstract. We prove global existence of non-negative weak solutions for a strongly 
coupled, fourth order degenerate parabolic system governing the motion of two 
thin fluid layers in a porous medium when capillarity is the sole driving mecha- 
nism. 



1. Introduction and the main result 

In this paper we study the following one-dimensional degenerate system of equa- 
tions 

fd t f= -d x [fd 3 x (Af + Bg)], 

[ d t g = -d x [gd* (/ + g)\ , 
which models the dynamics of two thin fluid threads in a porous medium in the 
absence of gravity. One of the fluids is located in the region bounded from below by 
the line y = and from above by the graph y = f(t,x), while the region occupied 
by the second fluid is located between the graphs y = f(t, x) and y = (/ + g)(t,x), 
f and g being non-negative functions. Furthermore, L is a positive real number and 
the positive constants A and B have the following physical meaning 

A , = fi+ Id + lw B . = 

M- Id /i-" 
We let /i_ [resp. //+] denote the viscosity of the fluid located below [resp. above], 
7uj is the surface tension coefficient at the interface y = f(t, x) between the wetting 
phases, while 7^ is the surface tension coefficient at the interface y = (/ + g)(t,x). 
The system (2.1) is supplemented by initial conditions 

/(0) = /o, 9(0) =90, x€(0,L), (1.2) 

whereby fo and go are assumed to be known, and we impose no-flux boundary 
conditions 

d x f = d x g = dlf = d 3 x g = 0, x = 0, L. (1.3) 
The system (1.1) has been obtain in [7], by passing to the small layer thickness 
in the Muskat problem studied in [6]. This is a widely used approach in the study 
of thin fluid threads because it reduces complex moving boundary value problems 
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to local problems defined by generally simpler equations. System (1.1) is strongly 
related to the Thin Film equation because if, for instance, / is constantly equal to 
zero, then g is a solution of the Thin Film equation 



when n = 1. We refer to the survey papers [1, 9] where many aspects concerning 
the Thin Film equation are discussed. It should be noted that similar methods to 
those in [7] have been used in [8] and [10] to rigorously show that, in the limit of 
thin fluid threads, the solutions of the moving boundary value problems for Stokes 
and Hele-Shaw flows converge towards the corresponding solutions (determined by 
the initial data) of the Thin Film equation (1.4), with n = 3 for Stokes and n = 1 for 
the Hele-Shaw flow. Compared with the Thin Film equation, system (1.1) is more 
involved because it is strongly coupled, both equations of (1.1) containing highest 
order derivatives of / and g, and, furthermore, there are two sources of degeneracy, 
because both / and g may be equal zero. Since both equations of (1.1) have fourth 
order, we cannot relay on maximum principles when studying problem (1.1). 
Corresponding to (1.1), we define the following energy functionals 



W,g):=\f \d x f\ 2 + -r^—\d x (f+9)\ 2 dx, £ 2 (f,g):= [ <D (/)+£<%) dx, 



whereby the function $ is given by 3>(s) := sln(s) — s + 1 for all s > 0. They will 
play the key role when constructing the weak solutions for the problem (1.1)-(1.3). 

Using these two functionals and Galerkin approximations, we prove that the prob- 
lem (1.1)-(1.3) possesses for non-negative initial data non-negative global weak so- 
lutions. To this end, we regularize first the system (1.1) and use the functional 
£2 to establish convergence of certain Galerkin approximations towards global weak 
solutions (of the regularized problem) which satisfy similar energy estimates as the 
classical solutions of (1.1)- (1.3). In a second step, we show that weak solutions of 
the regularized problem converge towards non-negative global weak solutions of the 
original system (1.1). The uniqueness of our weak solutions is left as an open prob- 
lem (this is still an open problem also for the Thin Film equation cf. [2, 11]). We 
note that it has been only recently shown in [13] (see also [3, 4]), in the context 
of the Thin Film equation, that the non-negative weak solutions found in [2] con- 
verge exponentially fast in H 1 towards flat equilibria. In our case, this is a further 
open question. The second order version of (1.1), when the fluids are driven only by 
gravity and surface tension is neglected, has been recently investigated in [5] where 
existence of non-negative global weak solutions which converge exponentially fast in 
L2 to flat equilibria is established (see also [7]). 

In order to state our main result, we introduce now the function spaces we work 
with. For each m £ N, we let H m := H m ((0, L)) be the L 2 — based Sobolev space 
and we let denote the closed subspace of H m which has {(pk ■ k £ N} as its 
basis. Herein, 



dtg + d x (g n d 3 x g) = 0, 



(1.4) 





00 := \A/-k an d 4>k '■= \/2/Lcos(k7rx/L), k > 1, 
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are the normalized eigenvectors of the operator —d 2 : H 2 — > L 2 with zero Neumann 
boundary conditions. To be more precise, / G if and only if the Fourier series 
associated to / converges towards / in H m . It is well-known that iJ^ = H 1 and, it 
is not difficult to see that, for m > 4, the boundary conditions (1.3) are satisfied by 
functions from this space. 

Given T G (0, 00], let Qt '■= (0, T) x (0,L). The main result of this paper is the 
following theorem. 

Theorem 1.1. Let fo,go 6 H 1 be two non-negative functions. There exist a global 
weak solution (f,g) of (1.1) with (/(0), <?(())) = (/oi5o) an d having the following 
properties: 

(1) / > and g > in (0, T) x (0, L), 

(2) f,g G L^iO^; H 1 ) 11 L 2 (0,T; Hi) 11 C([0,T],C a ([0, L])) for some arbitrary 
a G (0, 1/2) and Vfd 3 x (Af + Bg), ^gd 3 x (f + g) G L 2 (Q+), where 



Q+ :={(i,x) GQ T : (fg)(t,x)>0}, 



(3) 



/ f(T)tp dx - [ M dx+ f {Ad 2 J + Bd 2 x g)(d x fd x TP + /<9» dxdt = 0, 
io Jo JQ T 

L f(T)tjj dx - [ L M dx + f (d 2 x f + 8 2 x g) {d x gd x ^ + gd 2 x ^) dxdt = 







for all T > and -0 G H^. Furthermore, the weak solutions satisfy 

(4) ||/(r)|| Ll = |l/olk and \\g(T)\\ Ll = \\g \\ Ll , 

(5) £ 2 (f(T),g(T))+ [ (A-B)\d 2 x f\ 2 + B\d 2 x (f + g)\ 2 dxdt<£ 2 (f Q ,g ) 

Jq t 

for allT G (0, 00) , anc? 

(6) S 1 (f(T),g(T))+ [ f\d x (Af + Bg)\ 2 + Bg\d 3 x (f + g)\ 2 dxdt < £i(/ , 5 o) 



/or almost all T G (0, 00) 



We remark that since f(t) and g(£) belong to H 2 ^ for almost all t > 0, they satisfy 
homogeneous Neumann boundary conditions at x = and x = L for all such t. 

The outline of the paper is as follows: in Section 2 we introduce a regularized 
version of (1.1) and use Galerkin approximations to find, in the limit, global weak 
solutions of this regularized problem (see Proposition 2.1). Introducing the regular- 
ized system allows us on the one hand to use the energy functional E\ when dealing 
with the Galerkin approximations, and, on the other hand, to control the solutions 
of the regularized problem when they become negative. In Section 3 we show, by 
combining energy estimates for both functionals E\ and £2, that the weak solutions 
of the regularized problem converge towards non-negative global weak solutions of 
our original problem (1.1)-(1.3). 
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2. The regularized system 

In order to prove the Theorem 1.1 we shall regularize system (1.1) and use Galekin 
approximations to build global weak solutions for this regularized problem. These 
solutions are shown later on, in Section 3, to converge towards weak solutions of 

(1.1) . To this end, given e G (0,1], we define the Lipschitz continuous function 
a £ : K — > K by the relation 

n / s + e, s > 0, , v 

a£(s) := ( e> S < 0. ^ 

Furthermore, we define the convex function $ £ : K — )• K with 

{(s + e)ln(s + e)-(s + e) + l, s > 0, 
s 2 (2 2) 

— + sln(e)+eln(e)-e + l, s < 0. V ' 7 

Since we choose e < 1, it is easy to see that $ e (s) > for all s € R and that 
<3?" = l/a e . With this notation, we introduce the following regularized version of our 
original problem (1.1) 

-d x [a £ (fe)d 3 x (Afe + Bg £ )}, 

a { f \f&f f , M (t,x)G(0,oo xfi. (2.3) 

Of course, this system is coupled with the initial and boundary conditions (1.2) and 
(1.3). Compared to (1.1), the only difference is that we replaced at one place / and 
g in (1.1) by a e (f) an d a E (g), respectively, and penalize in this way the functions 
f e ,g £ when they take negative values (see the definition of a e ). Furthermore, by 
choosing the regularization in this way, we may still use the functional £\ to obtain 
useful estimates for the solutions of (2.3). For the problem consisting of (2.3) and 

(1.2) -(1.3) we prove the following result. 

Proposition 2.1. Let /o,<?o £ H 1 be two non-negative functions and e G (0,1]. 
There exist globally defined functions f e and g e with / e (0) = /o, <? E (0) = go and 
having the following properties: 

(i) Given T > 0, the functions 

fe, 9e G £co ( [0, T] , H 1 ) n L 2 (0, T; H 3 A ) n C{ [0, T] , C a ( [0, L] )) 

for some arbitrary a G (0, 1/2). 
(ii) For all T > and tp G H 1 we have 

-L r L 




f f e (T)^dx- f Mdx= [ a £ (f £ )d 3 x (Af £ + Bg £ )d x ^dxdt, 
Jo Jo Jq t 

/ 9e(T)^dx- g ^jdx= a £ {g £ )dl{f £ + g e )d x ip dxdt. 



o 
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(Hi) The following energy estimates are satisfied: 

pL pL pL pL 

(a) / f £ (T)dx= / f dx and / g £ (T)dx= \ g dx, 



JO JO 



(b) [ $ £ (f £ (T)) + B<5> £ (g £ (T))dx 
Jo 

+ [ {A-B)\dlf £ \ 2 + B\dl{f £ + g £ )\ 2 dxdt 
JQt 

< f <t> £ (f ) + B$ £ (g )dx 

J 

for all T £ [0, oo), and 
(c) £i(f £ (T),g £ (T)) 

^— [ a £ {f £ )\dl{Af £ + Bg £ )\ 2 + Ba £ {g £ )\dl{f £ + g £ )\ 2 dxdt 
~ B J Qt 

< £i(fo,go) for almost all T 6 (0, oo). 



+ A 



We will construct the global solutions of (2.3) by using Galerkin's method. In a 
first step we will find, by using the Picard-Lindelof theorem, Galerkin approxima- 
tions for the solutions of (2.3) which are defined on a positive time interval. Using 
the energy functional £±, we show then that in fact the approximations are defined 
globally. In a second step, we prove that the Galerkin approximation converge to- 
wards global solutions of the regularized system which satisfy energy inequalities for 
both energy functionals £\ and £2. Though /o and go are non- negative, it is not clear 
whether f £ and g £ preserve this property in time. However, we will show later on, in 
Section 3, that, for e — > 0, f £ and g £ converge uniformly to non-negative functions. 

2.1. Global existence of the Galerkin approximations. 

Given fo,go in H , the initial conditions of (1.1), we consider their expansions 

00 00 

/o = E fok^k, 9o = ^2 9ok< ^ k in 

k=0 k=0 
and, for each n € N, the partial sums 

n n 

fo := E ^ ofc ^ fc ' 9o := E 90k4>k- 
k=0 k=0 

We first seek continuously differentiable functions 

n n 

fe ■■= E F ^k, 9e ■■= E 

k=0 k=0 

which solve (2.3) when testing with functions from the vector space {<po<, • * • ? 4*n) 1 &fid 
additionally 

/ £ n (0) = /o n , 9e(0)=9o- 
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By construction, the functions (f™,g £ ) satisfy the boundary conditions (1.3) and, if 
we test (2.3) with constant functions, it follows at once that necessarily F £ and G £ 
are constant functions 

= /oo, G° £ (t)=g 00 , t>0. (2.4) 

Moreover, the tuple (F™, G%) := (F 1 , . . . , F™, G\, . . . , G%) is the solution of the initial 
value problem 

where ^ := (\I>i, * 2 ) : K 2n ->■ R 2n is given by 



{F n^ G ny = WjG n) 9 (i?,G?)(0) = (/oi, . . . ,/ n,501, • • • , <70n), (2.5) 

p2n , TTj2n 



xi4>i dl<j) k d x (t)j dx 



n pL / n 

^i,j(x,y) = ^2(Ax k + By k ) / a e I / 00 ^o + ^ 
fc=l ^° V !=1 

n r L ( n \ 

i &2,j(x, y) = Yl( Xk + Vk ^ J a£ ( foo^o + y^ 1 ) dl^kdx&j dx, 

for all x, y G M n . Since a £ is Lipschitz continuous, we deduce that $ is locally 
Lipschitz continuous on M 2n , and therefore problem (2.5) possesses a unique solution 
(F™,G™) defined on a maximal interval [0,T"). In order to prove that the solution 
is global, that is T£ = oo for all e G (0, 1] and n G N, we make use of the energy 
functional 6\. Indeed, since d x f £ ,d x g £ G (4>o, ■ ■ ■ , 4 > n)i we may use them as test 
functions for (2.3). Integrating by parts, we then get the following relation 

|ft</r.rf) 



1 



A-B 



1 



A-B 
1 

A-B 



L Ad x f?d t (d x f?) + Bd x f?d t (d x g n £ ) 

+ Bd x g n £ d t {d x fl l ) + Bd x g^d t (d x g^) dx 

Ad 2 J?d t f? + B [d 2 j:d t g n £ dx + d 2 x g n £ d t g n £ + %g?dtf?] dx 
L [Aa E (f?)%f?dl(Af? + B£) + Ba E (g?)d 3 J?d 3 x (f? + 9?) 



+ Ba e (f?)%£%(Af? + Bg n £ ) + BaM)&d x (f? + £)) dx, 

(2.6) 

and taking into account that £i(fQ,gn) < £i(fo,go) for all n G N, we find after 
integrating with respect to time that 



£l (f^ ng n {T)) + _L_ J a £ (f?)\dl{Af? + By 



n\\2 

Qt 

+ BaMMUe + 9s)\ 2 dxdt < £i(/ ,so) (2.7) 
for all T > 0. Whence, there exists a positive constant C, which is independent 
of time, such that |(F™(T), G"(T))| < C for all T < T e n . Together with (2.4), we 
conclude that for each n G N and e G (0,1], the Galerkin approximations (f £ ,g £ ) 
are defined globally. 
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2.2. Convergence of the Galerkin approximations. 

Let T > and e £ (0, 1] be fixed. From the energy estimate (2.7) we deduce that 

d x fe, 9x9e are bounded in Loo (0,T;L 2 ), (2.8) 



y/ajJf)dl{Af™ + Bg«), y/a^)^ + g n £ ) are bounded in L 2 (Q T ), (2.9) 

uniformly in n £ N and e £ (0, 1]. In view of a e > e and A > B, we obtain from 
(2.9) that 

d^f™, dig™ are bounded in L 2 (Qt), (2.10) 

uniformly in n £ N. Furthermore, by virtue of (2.4). we see that the mass of both 
fluids is preserved by the Galerkin approximations 

/ f£(t)dx= [ fodx and f g™(t)dx= f g dx for all t £ [0, T}. (2.11) 
Jo " J J •/ 

Invoking now (2.8), (2.11), and the Poincare-Wirtinger inequality we conclude 
that in fact 

g™ are bounded in L oo (0, T; H l ) uniformly in e £ (0, 1] and n £ N, (2.12) 

while, owing to (2.10) and (2.11), the same inequality implies 

g™ are bounded in L 2 (0, T;H 3 ) uniformly in n. (2-13) 

We consider now the partial derivatives with respect to time, and observe that the 
first equation of (2.3) can be written in the more compact form dtf™ = —d x H™ 
where, by (2.9), (2.12), and using the embedding H 1 L^, the right-hand side 
H™ := a{f^)(Ad 3 f r e l + Bd 3 g™) is bounded in L 2 (Qt) uniformly in e and n. Therefore, 
given £ £ H , we set 

n 
k=0 

and, using integration by parts, obtain 

\wm\o\ =\(d t mt)\c n )\ 
=m\d x Cn)\ 

<\\ H 7\\L 2 {QT)\\C l \\m 

<\\ H 7\\l 2 {Qt)\\^\\h^- 

This means that 

dtfe, d t g™ are bounded in L 2 (0, T; (H 1 )') uniformly in e and n. (2.14) 

Gathering (2.12)-(2.14), we obtain from Corollary 4 in [12], by making also use of 
the embeddings 

ff 1 ^' C^L})^^ 1 )' and H 3 ^ C 2+a ([0, L}) ^ (H 1 )' 
for a £ [0,1/2), that 

/ £ n , ^ are relatively compact in C([0, T], C a ([0, L])) n L 2 (0, T; C 2+a ([0, L])). 



8 B.-V. MATIOC 

Whence, for each e G (0, 1], there exist functions 



f £ , g £ G C{ [0, T] , C a ( [0, L) ) ) n L 2 (0, T; C^ +Q ( [0, L] ) ) 

and subsequences of (/") and (g £ ) (which we denote again by (/") and (g £ )) such 
that 

f?^f £ and g n £ ^g £ in C([0, T], C a ([0, L])) fl L 2 (0, T; C 2+a ([0, L])). (2.15) 

Moreover, we deduce from (2.13) that 

d p Je^dPf £ and - 8P x g E in L 2 (Q T ) for p = 1, 2, 3, (2.16) 

and therefore / £ ,cfe G L 2 ([0, T], i? 3 ). Additionally, since f™(t),g™(t) G H% we get, 
by virtue of (2.15), that d x f £ = d x g £ = at x = 0, L for almost all t G [0, T], which 
yields f E ,g E £L 2 ([0,T],Hl). 

2.3. Proof of Proposition 2.1. 

First of all, /™(0) = /q for all n G N and since f G H 1 we conclude that / e (0) = f 
for all e G (0, 1]. Similarly, we have g £ (0) = go f° r aJl e G (0, 1]. Furthermore, it is clear 
from (2.11) and (2.15) that the weak solutions (f £ ,g £ ) satisfy the relation (m)(a) of 
Proposition 2.1. 

We pass now to the limit in the energy estimate (2.7). By virtue of (2.9), (2.15), 
and (2.16) we have 

yfeU*j%{Aft + B&) - yfe(fJ%(Af e + Bg e ), 

m L 2 (Qt)- 

7^«)5 3 (/« + g n £ ) - v^)9 3 (/ £ + g £ ) 

Furthermore, by (2.15) we know that f £ (t) — > f £ {t) in H 1 for almost all t G [0, T], 
so that, by passing to the limit n — > oo in (2.7), we obtain the estimate (ra)(c) of 
Proposition 2.1. 

Claim (i) of Proposition 2.1 is now a simple consequence of the assertions (m)(a) 
and (ra)(c) of the same proposition. 

We now prove the assertion (ii) of Proposition 2.1. To this end, we pick an 
arbitrary function ip G H 1 and, testing (2.3) with ip n := X^o^l^)^' we obtain 
the following relations 

[ L f?(T)r dx - [ L f 0n r dx = [ a E (f?)d 3 x (Af? + Bg n e )d x r dxdt, 
Jo Jo Jq t 

L g n £ {T)r dx - [ L g 0n r dx= I aM)d 3 x (f? + g n £ )d x ^ dxdt. 



Jo Jq t 

Invoking (2.9) and (2.12), we see that a £ (f™)d 3 (Af™ + Bg£) and a £ (g™)d 3 (f™ + g%) 
are bounded in L 2 (Qt) uniformly in e and n. Using (2.15) and (2.16), we may even 
identify their weak limit 

a E (f?)d 3 (Af? + Bg £ ) - a £ (f £ )d 3 (Af £ + Bg £ ), 

in L 2 (Qt), (2.17) 

a £ {g n £ )dlU? + 9 n £ ) ^ a £ (g £ )d 3 (fe + 9 e) 
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and the assertion (ii) of Proposition 2.1 follows from the previous identities when 
letting n — > oo. 

We end this paragraph with the proof of the estimate (iii)(b) of Proposition 2.1. 
Let us observe that & £ (f™(t)) and $' e (f e (t)) belong to H 1 for almost all t <E [0,T], 
meaning that 



k <Pk, 



A:=0 



k=0 



(2.18) 

for almost all t E [0, T]. Of course, (2.18) is also valid when replacing / by g. In view 
of (2.18), we obtain the following relations 



dt 



[ L Mfe) + B^ £ {g n e ) dx = [ L *M)dtf? + B&M)dt9? dx 
Jo Jo 

n 

a £ (f?)d 3 x (Af? + Bfi) ^'eife)\4>k)d.. 



k=0 



+ Ba e (g n £ )dlU: + <£) £(**G£)l&)0*& dx, 

(2.19) 



k=0 



and, integrating with respect to time, we arrive at 
L <S> e (f?(T)) + B<S>M(T))dx 



a e {f:){Adlf: + B& x g») J2^e(f?)\fa)d x <t>k 

k=0 

n 

+ BaM)(d 3 Je n + dlg n e ) ^'MMk)d x <t>k 



k=0 



(2.20) 



dxdt 



+ f L $ E (fZ) + B<!> E (gZ)dx. 
Jo 



In order to pass to the limit n — > oo in relation (2.20) we have to determine what 
happens with the two integrals on the right-hand side of (2.20). Using (2.18), we 
have 

2 



Y^{®'sUs)\fa)d x <t>k-®"(fe)d*fe 



k=0 



L 2 {Qt) 



< 2 



x9k 



k=0 



L 2 (Qt) 



+ 2 



^(& E (f £ )\<j> k )d x cf> k -^(f E )d x fz 



fc=0 



MQt) 
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Taking into account that the first sum on the right-hand side of the latter inequality 
is the truncation of the Fourier series of ^(/ e n )<9 x /™ - $£(/ e )0 x /e, cf. (2.18), its 
norm may be estimated as follows 

2 



~ K(fe)\<f>k)d x 



k=0 



L 2 (Qt) 



■XQt) 

<2\\^(fe) - ^(f £ )\\l UQ J\d x f £ \\l 2 (Q T) 

+ 2|K(/?)HL 



(QtV 



\9xfl d x f £ \\ L2 
<2e 4 ||/ e n - / e |li 00 (Q T )||9 a ./ 6 || i2( Q T) 



(Qt) 



+ 2e 2 \\d x f £ - d x f E \\ L2 (Q T y 

We note that the last inequality has been obtained by using the fact that is 
Lipschitz continuous with Lipschitz constant e~ 2 and < < e" 1 , properties 
which readily follow from (2.1), (2.2), and the relation = l/a e . Invoking (2.15), 
we resume our calculation with 

2 

->n->oo 0. (2.21) 

k=0 



L 2 (Qt) 



£(*U/ e ")-^(/ e )l^0M* 

k=0 

Concerning the second term, we obtain from (2.18) that 

n 

^2(K(feMk)d x <P k -^(fe)d x f t 



k=0 



L 2 



k=n+l 



La 



for almost all i 6 [0, T], and Lebesgue's dominated convergence theorem yields 



Y,^'eU:)\4>k)d x( t>k-K{h)d x f l 

fc=0 

Gathering (2.21) and (2.22), we conclude that 

n 



l 2 (<9t) 



in L 2 (Qt) 



(2.22) 



(2.23) 



fc=0 



Clearly, (2.23) remains true if we replace / by g. We sum (2.17), (2.23), use (2.15) 
and the fact that both /o and g$ are non- negative to obtain from (2.20), when letting 
n — > oo, the desired assertion (iii)(b) of Proposition 2.1. 

3. The proof of Theorem 1.1 

We shall use the global weak solutions (f £ ,g £ ) of the regularized problem (2.3) to 
find, in the limit e — > 0, global weak solutions of our original system (1.1). The key 
role is now played by the second energy functional £2, which will be used to prove that 
the weak solutions we obtain are non- negative and to identify in -£2(0, T; H 2 ) a weak 
limit of the global solutions of (2.3). Using integration by parts, we may eliminate 
then from the right-hand side of (ii) Proposition 2.1 the third order derivatives of f £ 
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and g £ , for which we don't have any kind of uniform bounds, and obtain in the limit 
e — > the assertion (3) of Theorem 1.1. 

To do so, we collect first some estimates for the family (f £ ,g £ ) which have been 
already established in Section 2. We have to pay attention because some of the 
estimates proven before are uniform only with respect to n, and of no use in this 
final part. Invoking (2.9) and (2.12), we deduce the uniform boundedness of 

f £ , g £ in Lo^O, T;^ 1 ), (3.1) 

V^UDdKAf, + Bg e ), y^MjdUfe + 9s) m L 2 (Q T ), (3.2) 

while, by virtue of Proposition (m)(a) and (6), 

d x f £ , d x g £ are uniformly bounded in ^(Qr), (3-3) 

f f e (T)dx= [ fodx, [ g £ (T)dx= [ g dx, (3.4) 
Jo Jo Jo Jo 

for all T > 0. Lastly, we observe that the estimates (2.12) and (2.14) are both 
uniform with respect to e £ (0, 1] and n £ N. This implies that the families {/" : 
£ 6 (0, 1], n € N} and {g £ : e £ (0,1], n 6 N} are both relatively compact in 
C([0, T], C a ([0, L])), if a G [0, 1/2), and therefore 

(/ £ ), (g £ ) are relatively compact in C([0, T], C Q ([0, L])). (3.5) 

Consequently, there exist subsequences (f £k ) and (g £k ) and functions /, g such that 

f £k ^f and g £k ^g in C([0, T], C a ([0, L])), (3.6) 

while, owing to (3.1), (3.3), we conclude that f £ ,g £ are bounded in -£2(0, T;H 2 ), 
which ensures, after possibly extracting further subsequences, weak convergence in 
Ij2(Qt) °f the spatial derivatives up to order 2 

d p J £k ^dlf and dlg £k ^d p x g in L 2 (Q T ) for p = 1, 2. (3.7) 

Recalling Proposition 2.1 (i) and (3.6), we deduce that f,g £ L2(0,T; H^) for all 
T > 0. Moreover, the sequences (f £k ) and (<7e fe ) converge strongly towards / and g, 
respectively, in a different norm than in (3.6). 

Lemma 3.1. Given T > 0, we have: 

f £k ^f and g £k ^g in L 4 (0, T; H 1 ). (3.8) 

Proof. We prove only the assertion for /. Since f(t) and f £k (t) belong to for 
almost all t G [0, T], we conclude that their first order derivatives at and L must 
vanish. Whence, using integration by parts, we get 

£ (£\d x (fs k -f)\ 2 dx^j dt= £ (£d 2 x (f £k -f)(f £k -f)dx^j dt 

< f \\dl{fe k ~ f)\\l 2 \\fe k ~ /Hi dt < L\\f £k - f\\ 2 Loo{QT) \\d 2 x (fs k ~ /)lli 2 (Q T ), 
and, together with (3.3) and (3.6), we get the desired conclusion. □ 
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Particularly, we obtain from (3.8), that f £k {t) — > f(t) and g £k {t) — > g(t) in H 1 
for almost all t G [0,T], and together with the estimate (3.1) we conclude that 
f,g e L^iT-H 1 ). Furthermore, / e (0) = f and g £ (0) = g for all e G (0,1], so 
that (3.6) yields /(0) = /o and g(0) = go- The estimate (4) of Theorem 1.1 follows 
by combining (3.6), the assertion (m)(a) of Proposition 2.1, and Lemma 3.2 below. 

We use now the energy estimate (iii)(b) of Proposition 2.1, to establish the asser- 
tion (1) of our main result Theorem 1.1. 

Lemma 3.2. The functions f and g found above are non-negative. 

Proof. Assume that there exists (T, xq) G Qoq such that f(T, xo) < 0. Since by (3.6) 
f £k — > f in C{Q T ), we conclude that there exists a constant 5 > and ko € N with 
the property that f £k (T, x) < —5 for all x G [0, L\ with \x — xq\ < 5 and all k > k$. 
We then infer from (2.2) that 

<M/ £fc (T, x)) = /efc 2 ( J fc ' X) + f £k (T, x) ln(e k ) + e k ln(e fc ) - e k + 1 > ^ 

for all x and k as above. This contradicts the assertion (iii)(b) of Proposition 2.1. 
Clearly, the argument is true when replacing / by g, and this proves the claim. □ 

In order to deduce the energy estimate Theorem 1.1 (5), we recall (2.2) and notice 
that, for all k G N, we have <& £k (f £k ) > § £k (f £k ), where 

$ ( \ = \ {s + e k )\n{s + e k )-(s + e k ) + l, s > 0, 
* e * w- \ e k ln(e k )-e k + l, s < 0. 

Givent G [0, T], the sequence ($ £k (f £k (t))) is bounded in C([0, L]) and $ £fe (/ £fc (t)) — )■ 
$(/(t)) pointwise on [0, L]. Lebesgue's dominated convergence implies then 

hminf [ L <S> £k (f £k (T))dx > liminf [ $ Eh (f £k (T)) dx = f $(f(T))dx. (3.9) 
Jo fc ^°° Jo Jo 

Of course, the relation still remains true when replacing / by g. By virtue of (3.7), 
we may pass to liminf k ^oo in relation (iii)(b) of Proposition 2.1, and obtain in this 
way the desired energy estimate (5) of Theorem 1.1. 

To deal with the energy estimate (6) of Theorem 1.1, we observe first that for all 

k G N 

\a £k (f £k )-f\<e k + \f £k -f\, 

meaning, by (3.6), that 

a £k {fe k ) -> / and a £k (g £k ) -> g in C(Q T ). (3.10) 

For every positive integer m, we introduce now the set 

Qt '■= {(t,x) G Qt ■ f(t,x) > m~ l and g(t,x) > m" 1 }, 

where we may control, by virtue of the estimate (iii)(c) of Proposition 2.1 and (3.10), 
the third order derivatives of both f £k and g £k : 

(dlf £k ), {dlg £k ) are uniformly bounded in L 2 {Qt). (3.11) 
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Taking into account that = U m Q™ , we may assume, after possibly extracting a 
further subsequence, that 

dlf £k ^dlf, dlg Ek ^dlg mL 2 (Q?) 

for all m £ N, which, together with (3.10), implies 



a £k (f £k )d 3 x (Af £k + Bg £k ) - y^d^Af + Bg), 

in Li(Q^). (3.12) 

a £k (g £k )d 3 x (f £k +g £k )^ sfgdltf + g) 

In fact, by virtue of (3.2), the weak convergence in (3.12) takes place in L 2 (Q™). 
Recalling Proposition 2.1 (m)(c) and Lemma 3.1, for k — > oo, we obtain the desired 
estimates (2) and (6) of Theorem 1.1. 

In order to complete the proof of Theorem 1.1, we are left to prove the relations 
(3). To this end, we pick ift £ H\. Since a £k is Lipschitz continuous, we obtain from 
Proposition 2.1 (i) that a £k (f £k (t)) G H 1 for almost all t £ (0,T) and 

dx{a £k {f £k )){t,x) = X(o,oc)(fe k )d x f £k , a.e. in Q T , 

whereby X(o,oo) denotes the characteristic function of the interval (0, oo). Integrating 
by parts in the first relation of Proposition 2.1 (ii) , we arrive at 

[ L f £k {T)i,dx- [ L Mdx= [ T [ L a £k {f £k )di{Af £k +Bg £k )d x ^dxdt = h, k +I 2 , k 
Jo J J J 

(3.13) 

where 

h,k ■=- d x (a £k (f £k ))d 2 x (Af £k + Bg £k )d x i>dxdt, 
Jq t 

h,k ■=- a £k (f £k )dl(Af £k + Bg £k )d 2 x i) dxdt. 
JQt 

We note the use of ip G to eliminate the boundary terms in (3.13) due to 
dMO) = d x ^{L) = 0. 

Combining (3.7) and (3.10), we obtain for k — > oo that 

h,h -)• - / f(Ad 2 x f + Bd 2 x g)d 2 x ^ dxdt. (3.14) 
JQt 

We consider now the integral I\ and notice that in order to show the relation 

h,k ->• - / d x f(Ad 2 J + Bd 2 x g)8 x ^ dxdt (3.15) 
JQt 

it suffices to prove that 

d x (a £k (f £k )) ^ 8 x f in L 2 (Q T ). (3.16) 

To this end, we write d x (a £k (f £k )) - d x f = (d x (a £k (f £k )) - d x f £k ) + (d x f £k - d x f) , 
and conclude from Lemma 3.1 that {d x f £k — d x f) — > in L 2 {Qt)- Furthermore, the 
first term may be written as {d x (a £k {f £k )) - d x f £k ) = (x(o,oo)(/e fc ) ~ l) d x f £k , and 
since d x f £k — > d x f in L 2 (Qt), there exists a function F € L 2 (Qt) such that, after 
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possibly extracting a further subsequence, \d x f £ .\ < F almost everywhere in Qt (see 
the proof of Theorem 3.11 in [14]). We show now that (x(o,oc)(fs k ) ~ l) dxfe k ~^ 
almost everywhere in Qt- Indeed, since d x f ek —> d x f in L,2(Qt), we deduce that 
9xfe k ~~ ^ almost everywhere on the set [/ = 0]. Furthermore, on the set [/ > 
0], relation (3.6) implies pointwise convergence (x(o,oo)(/efc) ~~ -U ~ ^ 0. Lebesgue's 
dominate convergence theorem implies now the desired relation (3.16), and implicitly 
(3.15). 

To conclude, we sum (3.6), (3.14), (3.15) and let k — > oo in relation (3.13) to 
obtain the first identity of Theorem 1.1 (3). The corresponding relation for g follows 
similarly. 
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